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TECHNICAL MEMORANDUM 1381 


ON THE DETERMINATION OF CERTAIN BASIC TYPES OF 
SUPERSONIC FLOW FIELDS* 


By Carlo Ferrari 
SUMMARY 


A quite universal mode of attack on problems which arise in super- 
sonic flow, whether connected with flow over wings or over bodies of 
revolution, is explained, first, in great generality, and then in more 
detail, as specific applications to concrete cases are illustrated. The 
method depends on the use of Fourier series in the formal definition of 
the potential governing the flow and in the getting up of the boundary 
conditions. This new formulation of the many problems met in supersonic 
flow is really an extension of the doublet type of "fundamental solution" 
to higher order types of singularity. The limitations and, in contrast, 
the wide field of applicability of such a means of handling these prob- 
lems with complex boundary conditions is discussed in some detail, and 
a specific example of a wing-body interference problem is cited as proof 
of the versatility of the method, because the results obtained by applying 
the techniques expounded herein agree well with experimentally determined 
data, even for the quite complex configuretion used to exemplify the kind 
of problem amenable to such treatment. 


1. INTRODUCTION 


For purposes of analytic treatment of the flow problem to be con- 
sidered here the usual rectangular Cartesian coordinate system is employed 
with the x-axis taken to lie in the direction of and having the same sense 


as the uniform (undisturbed) free-stream velocity, Vo. This free-stream 


velocity, Vo, is taken to be supersonic in the discussion that follows; 


i.e., Va> Co where C, denotes the velocity of sound in the undis- 
turbed stream. The flow of the gaseous fluid to be investigated is to 


¥ 
"Sulla determinazione di alcuni tipi di campi di corrente iper- 
sonora," Rendiconti dell'Accademia Nazionale dei Lincei, Classe di 
Scienze fisiche, matematiche e naturali, serie VIII, vol. VII, no. 6; 
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be considered as resulting from the superposition upon the free-stream 
velocity of a nonuniform flow, having velocity components that are des- 
ignated as Vy, Vy, and Vz, and lying in the direction of the respec- 
tive axes (x,y,z) of the coordinate system. This nonuniform super- 
imposed flow is supposed to be small enough, in comparison with the 
speed of sound, Ce, that it is permissible to neglect the ratios Vy/Cw; 
Vy/Coos etc. in the equations governing the flow. 


It is taken for granted that, under the conditions stated above, 
there exists a velocity potential describing the flow in question, and 
in practically all cases which are of any interest for actual designs 
it will really be true that this assumption can be made legitimately. 
If it is then agreed that the nonuniform superimposed part of the flow 
is to be denoted by the potential Ø, it will be recognized that this 
potential will have to satisfy the relationship: 


a - u2) oe SE + of. x (1) 


where the free-stream Mach number, M,, is defined as M, = Vof Con and, 


of course, the potential Ø will also have to obey the boundary condi- ` 
tions which are peculiar to each stated problem. 


A way of handling the determination of the function Ø, so that it 
will satisfy equation (1) and so that it will obey the imposed boundary 
conditions, will now be explained, and its usefulness illustrated by 
consideration of problems which can be attacked by this means, both in 
the case of lifting surfaces (that is, wings) as well as in the case of 
bodies of revolution. The proposed method is based on the use of Fourier 
series. Although this technique does not afford complete universality 
in treatment of all the posed problems, as will be more clearly pointed 
out in what follows, it can be used to fine advantage in a goodly number 
of situations by replacing the procedures which are based on the Fourier 
or the Laplace transforms (which, for that matter, have just as restricted 
limits of applicability as the analogous ones which arise in connection 
with the approach being discussed herein) or by being substituted in 
place of the techniques which stem from use of the "fundamental" (source, 
sink, doublet) solutions to equation (1), or from use of transformations 
carried out in the ‘complex plane. 


2. PROBLEMS HAVING TO DO WITH FLOW OVER WINGS 


As usual, the wings are imagined to be very slender and so placed a 
that the wing span lies along the y-axis; i.e., the long dimension is 
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out the y-axis (see fig. 1). Let the equations which define the ventral 
and dorsal surfaces of the wing surface be, in fact, given in the form: 


z= 2. (x,y) and z = z,(x,y) 


and then the slenderness of the wing is supposed to be slight enough 
that the above-defined values of z will be so small at all locations 
on these surfaces as to make it possible to accept the fact that the 
derivative gz is, to all intents and purposes, equal to the direction 
cosine, with respect to the free-stream x-axis, of the normal to the 
surface. It is further assumed that the wing is immersed in a stream 

of supersonic flow which has a constant value for its component lying 
in the direction of the x-axis, of magnitude Vo. The component in the 
direction of the z-axis, meanwhile, is assumed to be known, but of rela- 
tively small size in comparison with the V, velocity, and it may take 
on various values, which will be denoted by V,'. If, now, the potential 


describing the flow perturbed by the wing is denoted by d this poten- 
tial will have to satisfy equation (1), and it will also have to conform 
to the conditions which are imposed at the boundaries. These further 
(boundary) conditions may be stated as follows: 


(2) Upstream of a certain surface, which may be immediately defined 
Just as soon as the wing-like body is specified which is to invest the 
impinging stream, the value of g is zero; i.e., the basic condition is 


d=0 (2) 
(3) On the wing surface, it must be true that 


($) = We E ~ vy" cos (8,2) = E(xy) (3) 


wherein the value of z to be employed is either the zy or za quanti- 


ties, depending on whether one is concerned with a point which is lying 
on the under ventral surface or on the upper dorsal surface, respectively. 


The notation cos (2,2) signifies the cosine of the angle between the 
z-axis and the unit vector taken in the direction of the exterior normal 
to the wing surface in question; i.e., this vector is, represented by the 
vector n, and under the present hypothesis cos (2,2) = tl. 


4 NACA. TM 1381 


It is convenient to distinguish between two basic types of problem 
which come under this kind of analysis, and to make the differentiation 
on the basis of the sort of boundary conditions met with in each type; 
that is, 

symmetric Types of Configuration 


In this case, the boundary conditions to be satisfied on the wing 
may be expressed in the form 


(22) = H(x,y) (3') 
z=0* 


dz 


and 


dp _ | 
Ea = -H(x,y) 


Asymmetric Types_of Configuration 
In this case, the boundary conditions are expressed as 


The first type of problem corresponds to a configuration for which 
the wing has a zero angle of attack with respect to the free-stream 


undisturbed flow, Vo, and which possesses a symmetric profile. The sec- 
ond type of problem corresponds to a configuration for which the wing 
is a flat plate, but which has any local angle of attack whatsoever, 


with respect to the free-stream vector, Vo, so Long as it is small. 
5. DEVELOPMENT OF THE CASE OF THE SYMMETRIC TYPE OF CONFIGURATION 


In this case it will suffice to examine the flow solely in the 


upper half-plane, where z >0. If GI) (2,72) stands for the flow 


which takes place in this upper region, anà if b(2) (x,y,z) represents 
the flow in the nether région, then, of course, 


6°?) (x,y,2) = g(t) (0,5, -2) 
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The boundary conditions in this case are composed of equations (3'), 
together with the restriction that 


(1) 
pe = 0 (for locations lying beyond the region 
02 / 7-0 occupied by the wing surface) (2') 


Now let the definition of the function describing the velocity com- 
ponent at the wing surface, and also the potential function itself, be 
cast into the convenient forms 


E*X(x,y) = Vo > Hm (8) cos emy = Vo S ml) cos 5 IM 
Ti m 


and 


GL) (x,y,z) = g(t) = Vo b Do d_ (8,6) SOS = mq (4) 


wherein & = x/b, n= y/b, and € = z/b, while b is a suitable length 
used for purposes of nondimensionalization. The value used for b will 
be equal to the semispan of the wing in the case where the leading edge 
of the wing is supersonic everywhere, and provided that the wing tips 
are cut off in such a way that the wing surface remains outside of the 
tip Mach cones emanating from either one of the wing-tip extremities out 
at the farthest reaches of the wing span. The value used for b will 
be larger than this semispan Just defined, if, in contrast, these geo- 
metrical relationships do not hold; the magnitude employed for b in 
this latter case is illustrated in figure 2. 


Finally, it should be observed that H* is a periodic function 
of y, which is equal to the values taken on by the function H at the 
wing's surface and it is zero for points lying out of this region, and 
this definition is to hold throughout the spanwise interval for which 
-b Sy b. 


The fact that it is possible to write H*(x,y) in the form given 
as equation (4) (i.e., the possibility of expressing the component- 
velocity field describing the normal velocities to the wing surface by 
means of a Fourier series instead of in terms of a Fourier integral) 
stems from the property already noted to the effect that the perturba- 
tions, which are created at any arbitrary point P(x,y) whatsoever, do 
not make themselves felt anywhere outside of the Mach cone emanating 
from P. As a result of this situation, therefore, as far as the 


6 NACA TM 1381 


determination of the fiela of flow about the given wing is concerned, it 
makes no difference to this flow whether one considers the wing to be 
operating by itself as an isolated entity within the impinging stream or 
whether, instead, one imagines it to be accompanied by an infinite number 
of reflections of this primary wing in the planes y = tmb. 


If one now inserts the second of the expressions given as equa- 


tion (4) into equation (1), it will be seen that this differential equa- 
tion reduces to 


fm _ no dm _ 2 
le «Brn (5) 


wherein B2 = Mf - 1 and wheré k replaces the constant a 


Meanwhile, it is also evident that, on the basis of the first of 
the formal developments given as equation (4), the boundary condition 


reduces to 
oe 3 


The expression given as equation (5) above is formally analogous to 
the so-called "telegraph equation,” and its solution, which is suitable 
for applying the type of boundary condition ‘exemplified in equation (6), 
is 


E -BÉ 
Bn = M bm(6') Jol V(E-81)2 - Ta de (7) 


where Jog is the cylindrical Bessel function of zeroth order. 


Consequently, the vertical derivative turns out to be 


od EEG k 2 o i FE 
= = xh (£-Bt) + rkt (6°) 31$ le-8')* - BC? | 
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and, because of the boundary condition. (6), it follows that 


so that the sought potential must have the form 
i E -B6 
fa = -F J. Em(e')Jo de' (7!) 


4. DEVELOPMENT OF THE CASE OF THE ASYMMETRIC TYPE OF CONFIGURATION 


+ 


The possibility of being able to find solutions to such asymmetric 
problems by means of the method being propounded here is restricted in 
this case to those configurations for which the leading edge as well as 
the trailing edge of the wing are supersonic, and where the wing tips 
are cut off in such a way that the wing surface lies outside of the Mach 


cone emanating from the very tip of the leading edge where the maximum 
span occurs. 


Under these circumstances the boundary conditions are constituted 
from the restrictions given as equations (3"), and of equation (2') 
once again. If one then follows the same procedure as was utilized in 
section 3, it follows that the expression for the sought potential is 
formally given as (ref. 1) 


¢-B|t| È a ag! 
pe _ k ‘ya (Kite -e2 — peH 
Dn rb | Bel + Ù cf até "BE = 8") W(e-e')?-Be¢? 


(8) 


where hm is, a priori, an undetermined function, and where it should 


be recognized that the + sign is to be employed for the lower half-plane 
where € < O, and where the - sign is to be employed for the upper half- 
Plane where € > 0. It is evident, therefore, that the derivative of g 


with respect to € will be continuous along the plane € = 0, but the 
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derivative of @, with respect to È will be discontinuous, and the 
"jump" will be of such size that 


(se t=07 (2 t=07 


It is clearly permissible here again to concentrate attention solely 
upon the disturbed flow in the upper half-plane where È > O, therefore, 
because the observation just made above will tell one how to compute what 
the flow will be in the other lower half-plane, once the former is 
obtained. 


holds true. 


The boundary conditions in this instance may now be recast into the 
form 


. É 1 
xBhy(£ ) + / hy(&') Jo (x Et) de! - 
P ia 
af bm(6') Jz ( bi E = Gylé) 


provided, as in the previous section, one sets up the convenient con- 
vention that G*(x,y) is to represent a periodic function in y that 
is to be equal to the values taken on by the function G(x,y) at the 
wing's surface, and it is to be zero for points lying out of this region. 
This definition is to hold throughout the spanwise interval for which 

-b Sy Sb. In addition, the form of G*(x,y) is to be assumed, specif- 
ically, to have the appearance 


(2,7) = Vo D. Glè) cos fm 


while it has also been assumed that the derivative of a function by the 
sole parameter upon which it depends is to be denoted by a dot over the 
function, that is, 


Liga ® 
Palë) = & 


4 
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The integro-differential equation defining h, may also be immedi- 
ately simplified to the compressed expression 


. 5 
sBhm(£ ) +85 bm(8') (Jo + Jo) di' = Gp(£) (9) 


Now apply a Laplace transformation to this integro-differential 


equation (i.e. , multiply through by the factor e` PS and integrate 
from O to co) . Thus, one obtains 


where a bar over a symbol serves to indicate that this quantity stands 
for the Laplace transform of the function so designated. 


Standard tables of Laplace transforms could be consulted to check 
these results, which may now be simplified by noting that 


2 
ke 2 2 2 2, ke 
+ = - + + + - + 
Pte) atte? ft 


t = 
sa ieee Ke fre + 
Be p2 Be Be Be 


2 


l 
MD 
to td 
AN 
td 
ho 
+ 
Wol h 
I 
nr La 


Thus the Laplace transforms of equation (9) simplifies to 


cook + ahha + - ) hm = Gn 


or the explicit expression for tbe Laplace transform of the unknown func- 


tion hy, is given in the form 
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so that finally one may invert the transformation to obtain 


SEE nle o [È (e - 8) de! (10) 


Once having obtained the value of hm, it is easy to write down the 


expression for the component of velocity lying in the x-direction and 
located at the wing-surface, because one has simply that this component 

is given by the partial derivative of the potential Ø, taken with respect 
to €, and evaluated at the plane of the wing; i.e., one finds that 


A) _ ZA À 
ASI a oe ee 


Furthermore, the formula giving the lift on the wing is just 
ml 
2,2 5 lp(2/0) 
L = 8p,V_b ) (1) = —— (11) 


for m= 1,5;5,. » a, etc. 


where the symbol, 2, is used to denote the distance along the x-axis 
measured back from the leading edge of the root chord to the projection 
into the plane of symmetry of the trailing edge of the wing-tip profile. 


In regard to the moment taken about the y-axis, it is apparent that 
it may be computed from the relation: 


l/b 
My Erp DI > a COB fom 3) T En (ë )ag 
m - 


m+ 


2rp bv 2 È (E (1) ° [Ehm a f° balé) del (12) 


for m = 1,3,5,. . ., etc. 
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5. PROBLEMS HAVING TO DO WITH FLOW PAST BODIES OF REVOLUTION 


The procedure discussed in the preceding sections can be extended at 
once to apply also to the solution of problems which are concerned with 
the flow over bodies of revolution. 


For this purpose let a cylindrical coordinate system (x,Y,9) be 
set up, and then the equation which governs the potential, Ø, being 
sought will have the form 


6 LL, aF pd 
52 * ¥ ox * He * s2 (13) 


Now let R = R(x) be the equation of the meridian line of the body, 
and let it be assumed that R is sufficiently small at all locations 
along the body so that the direction cosine of the normal to this merid- 
ian line, measured from the Y-axis, may be taken to be equal to unity. 


Furthermore, let Vne! stand for the component, taken in a direc- 


tion perpendicular to the circular cross-section of the body, which 
arises from the impinging flow which invests the body. Then the boundary 
condition which must be satisfied at the surface of the body may be 
expressed mathematically by the relation 


For sake of simplicity, it is also now assumed that the treatment 
to be developed is to be restricted to the case where symmetry with 


respect to the semiplanes 8 = t5 exists in the incident flow. Under 


this hypothesis it is convenient to write the normal velocity components 
and the potential being sought in the following explicit formulations: 


We eie) T (©) sin m 9 
m 


(14) 


oe. 
Il 


Vo © Fn Y)Y" sin m 6 
m 
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If one now inserts the second of the expressions given as equa- 
tion (14) into the differential equation governing the flow (13), it 
will be found that the defining equation for the potential will have 
the form 


. 
MELE i sy =B L (15) 


and the boundary condition turns out to be 


3 (n) = -Fn (15') 


oY flyer 


A suitable sdlution to equation (15), which can be made to satisfy 
the boundary condition being imposed as equation (15'), will be found 
to be 


tn = (2 2)" % (16) 


where 


O 
do = J È fo(x-BY cosh u)du 
arc cosh = 


= 
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Thus, the successive individual potentials are given by the expressionst 


O 
Yø = -B 1 a f1(x-BY cosh u) cosh u du 
arc cosh SY 
O 
Yio = BÊ È m fo(x-BY cosh u) cosh® u du (17) 
arc cosh ay 
etc. 


Upon imposition of the requirement that the boundary condition (15') 
is to be satisfied, one obtains a set of integral equations which serve 


itranslator's note: It was pointed out on page 630 of an article 
by R. H. Cramer in the Journal of the Aeronautical Sciences, vol. 18, 
no. 9, September 1951, entitled “Interference Between Wing and Body at 
Supersonic Speeds - Theoretical and Experimental Determination of 
Pressures on the Body," that the result given here for p, for m> 1l, 
is incorrect; the correct formula is, for m = 2, 


O 
YG, = BÈ f fo(x-BY cosh u) (2 cosh® u - 1) du 
Vere cosh ay 


while, in general, the use of hyperbolic functions of muitiples of the 
argument u gives a more compact form, which is easy to work with; i.e., 
in general it is true that 


Yén = BH) T fp (x-BY cosh u) (cosh m u) du 


arc co = 
cosh =F 
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to determine the arbitrary functions fm, which are a priori unknown. 
Thus, applying these conditions, one finds that? 


O 
part (aye f fm(x-BR cosh u) (cosh uy@1 du = -Fn (18) 
arc cosh oA 


The determination of the values of the fn's appearing in formula (18) 


may be carried out by using a step-by-step procedure which is entirely 
analogous to the one employed by Von Kérmén in his work on determining 
the flow about a body of revolution at zero angle of attack. 


It is important to point out that if one only has in mind to caleu- 
late the force distribution along the axis of the body and the corre- 
sponding moment, and if one is not interested in knowing the local veloci- 
ties or pressures around the body, then it is merely necessary to 
calculate fọ and fp. 


6. APPLICATIONS 


The procedure that has been propounded above has been applied (ref. 2) 
to the situation arising in the study of the question of wing-body inter- 
ference. The wing-body configuration considered in this particular appli- 
cation of the method is depicted in the appended figure 3. 


The wing used in this configuration is a flat plate, whose semispan 
is equal to Ro, where Rọ is the radius of the circular cross-section 


taken through the body at the location where the body is widest. The 
leading edge of this wing is located 5Ro downstream from the tip of the 


nose of the body. The free-stream flow is impinging on the body at a 
speed which is twice the speed of sound in the undisturbed stream. 


“Translator's note: In view of the correction pointed out in Note 1 
above, it will be seen that this formula for determining the f, functions 
is also incorrect, except for m= l; for higher integral values of m, 
the correct formula is: 


ala h a. x Ém(x-BR cosh u) [cosh m u cosh u| du = -Fy 


arc cosh DE 


a 
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The curves shown in figure 3 give the value of the pressure coeffi- 


cient, acta, out along the span of the wing, in the mid-chord loca- 
Sora 

tion (i.e., along the wing axis), for points on the upper (dorsal) side 

of the wing. These coefficients have been calculated by the method out- 

lined in section l, and there are shown results for various angles of 

attack, which apply to such points on the upper side of the wing pro- 

files at their mid-chord positions. 


In addition, some experimental test points obtained by R. H. Cramer 
(see ref. 2) are also plotted on these curves. These results were 
obtained from experiments carried out in the supersonic tunnel of the 
Daingerfield Aeronautical Laboratory. 


The agreement between the computed and experimentally determined 
results is very good from a qualitative viewpoint. In regard to the 
more precise details of the quantitative comparison between the results 
it is worthy of note that the experimental results exhibit a certain 
amount of dissymmetry as one passes from positive angles of attack to 
negative angles of attack. Such a dissymmetry cannot be predicted, or 
should not be expected, from the type of theoretical treatment being 
considered here. 


In order to bring about a more valid comparison of these results, 
it would appear logical, in face of such evident dissymmetry, to take 
for the representative experimental value, at a given value of the angle 
of attack, B, the one which is obtained from averaging the result 
obtained for an angle of attack equal to +8 with the result obtained 
at -B. Such average values have been computed and are designated in 
the plots of figure 3 by means of solid circles. These adjusted values 
lie much closer to the theoretically derived curves at almost all 
locations. 


+ 


Translated by R. H. Cramer 
Cornell Aeronautical Laboratory, Inc. 
Buffalo, New York 
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Figure 1.- Orientation of coordinate axes and location of typical wing 
plan form therein. 


Figure 2.- Definition of the interval of periodicity required for application 
of the Fourier series technique when leading edges are subsonic. 
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Figure 3.- Pressure distribution along the wing axis: Comparison of 
experimental results with predictions based on the method expounded 
in section 4. 
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